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Abstract 

We present the spatial distributions of dust and polycyclic aromatic hydrocar- 
bons (PAHs) in the eUiptical galaxy NGC 4125, revealed by AKARI and Spitzer. 
NGC 4125 is relatively bright in the dust and the PAH emision for elliptical galaxies, 
although it certainly possesses diffuse interstellar hot plasma indicated by the high 
spatial resolution X-ray data of Chandra. We investigate how the dust and PAHs 
interact with the X-ray plasma or avoid the interaction by comparing their spatial 
distributions. We find that the distributions of the PAHs and dust are different from 
each other, both showing a significant deviation from a smooth stellar distribution. 
The PAH emission predominantly comes from a dust lane, a compact dense molecular 
gas region in the galactic center, where the PAHs are likely to have been protected 
from the interaction with the X-ray plasma. The dust emission has more extended 
structures similar to the distribution of the X-ray plasma, suggesting their interaction 
to some extent. We also discuss a possible origin of the dust and PAHs in the galaxy. 

Key words: dust, extinction — infrared: galaxies — galaxies: individual 
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1. Introduction 



Many elliptical galaxies possess an appreciable amount of X-ray-emitting hot plasma 
providing a harsh interstellar environment for the survival of polycyclic aromatic hydrocarbons 
(PAHs) and dust grains. Dust with a size of 0.1 /im is destroyed through sputtering on such 
a short timescale as 10^ yr for the plasma density and temperature of 10~^ cm~^ and 10^ K 
typical of elliptical galaxies (Draine & Salpeter 1979; Tielens et al. 1994). Recently, Micelotta 
et al. (2010) showed that the destruction timescale of typical PAHs (50 C atoms, size~ 6A) m 
hot plasma with temperature higher than 3x10^ K is so short as 10^ yr for the above plasma 
parameters, which is 10^~^ times shorter than lifetimes for dust with the same size. Thus the 
presence of dust grains in X-ray elliptical galaxies would imply that they were supplied very 
recently and the presence of PAHs would even suggest that there are some ways for them to 
avoid interaction with the hot plasma. 

Despite such a hostile environment, it has been found that a significant fraction of X-ray- 
emitting elliptical galaxies contain a considerable amount of dust, which cannot be explained 
solely from replenishment by old stars (e.g. Knapp et al. 1989; Goudfrooij & de Jong 1995; 
Temi et al. 2004; Temi et al. 2007). Surprizingly some of them even show the presence of PAHs 
(e.g. Kaneda et al. 2005, 2008a). Any clear signature of interaction of the dust and PAHs with 
X-ray plasma has not yet been obtained either from the relations of their emission luminosities 
among galaxies or from their spatial distributions within a galaxy. The former probably implies 
that the dust and PAHs are supplied by transient events rather than continuous events, their 
amounts depending on how recently the events happened. The latter is mainly due to a paucity 
of spatial information on the dust and PAH emission. For NGC 4589, the deep mid-IR spectral 
mapping with the Spitzer/IRS has revealed the distinctive spatial distribution of the PAH 
11.3 /im emission perfectly following the minor-axis dust lane (Kaneda et al. 2010). The 
distributions of the far-IR dust continuum and the PAH 17 /im emission are extended more 
widely in directions similar to each other and entirely different from the stellar distribution. 
Although the galaxy is detected in the X-ray, it is so faint that the distribution of diffuse hot 
plasma to be compared with those of the dust and PAHs cannot be derived. 

In this paper, we report on the mid-IR spectral mapping and near- to far-IR imaging 
observations of NGC 4125, an X-ray elliptical galaxy categorized as E6pec (Filho et al. 2002). 
The mid-IR spectral mapping data are obtained with the Infrared Spectrograph (IRS; Houck 
et al. 2004) onboard Spitzer (Werner et al. 2004), while the imaging data are obtained with 
the Infrared Camera (IRC: Onaka et al. 2007) and the Far-Infrared Surveyor (FIS; Kawada 
et al. 2007) onboard AKARI (Murakami et al. 2007). NGC 4125 has an intermediate X-ray 
luminosity (Log Lx = 40.94 erg s~^) among the ROSAT galaxies cataloged by O'SuUivan et al. 
(2001). The galaxy habors an emission line nucleus with a low-ionization nuclear emission- line 
region (LINER) spectrum (Heckman, Balick, & Crane 1980; Willner et al. 1985). Schweizer & 
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Seitzer (1992) concluded that NGC 4125 is a (6 — 8) x 10^ yr old remnant of a major disk-disk 
merger based on the UBV colors and the amount of merger-induced fine structures. Indeed 
the galaxy shows ionized gas emission along the major axis, which does not follow simple 
circular motions (Bertola et al. 1984). The kinematic profiles of stars suggest that there is 
a colder substructure in the innermost center (Pu et al. 2010). Verdoes Kleijn & de Zeeuw 
(2005) found a dust lane of size lf'3 in the center of the galaxy. Contrary to the suggested 
link between merging and X-ray brightness, the galaxy is not so X-ray-bright as hot gas-rich 
elliptical galaxies (Fabbiano & Schweizer 1995). Yet NGC 4125 certainly possesses diffuse 
interstellar hot plasma as indicated by the high spatial resolution X-ray data of Chandra and 
XMM-Newton (Humphrey et al. 2006; Fukazawa et al. 2006) and an appreciable amount of 
far-infrared (IR) dust (Knapp et al. 1989; Goudfrooij & de Jong 1995; Temi et al. 2004; Temi 
et al. 2007), thus is an ideal elliptical galaxy to search for signature of interaction of the dust 
and PAHs with the X-ray plasma. 

The contents of the paper are based on our new datasets of the IRS spectral maps 
and the AKARI images. NGC 4125 was also observed within the framework of the Spitzer 
SINGS legacy program (Kennicutt et al. 2003), although details of the individual galaxy have 
not been reported in a paper. As for the IRS observations, we performed spectral mapping 
observations, which were 2 — 8 times deeper and 2 — 3 times wider than the SINGS observations, 
depending on the spectrometer modules. Our AKARI image data provide new information at 
the photometric bands that the IRAC and MIPS do not have, especially at wavelengths of 11 
/im, 15 /im, and 90 /im. Below we adopt 22.2 Mpc for the distance to NGC 4125 (Tonry et al. 
2001; Jensen et al. 2003), where 1 kpc corresponds to the angular size of 9f'3. 

2. Observations 

The Spitzer/IRS spectral mapping observations were performed in May 2009, while the 
AKARI/IRC and FIS imaging observations were carried out in November 2006 and May 2007. 
The observation log is listed in table 1. The Spitzer data were taken in part of our Guest 
Observers (G05) program (PI: H. K.; program ID: 50369). Figure la shows the slit positions 
of the Short-Low (SL; 5.2-14.5 /im) and Long-Low (LL; 14-36 yum) modules. For ~ 1' x 1' area 
around the galactic center, we made a 33 x 2 grid for SL with a ramp duration of 60 s for each 
exposure, and a 13 x 2 grid for LL with a ramp duration of 30 s for each. We shifted the slit 
in directions perpendicular and parallel to its length by lf'85 and 3f'0, respectively, for SL, and 
by 5f'l and 9f'0, respectively, for LL. We repeated the same spectral mapping observations two 
times for LL. Starting with the basic calibrated data (pipeline version S18.7.0) reduced by the 
Spitzer Science Center, we extracted spectral image data by the CUBISM software (version 1.7; 
Smith et al. 2007a) following standard procedures. From the spectral image data, we created 
the spectra of the central 15" x 15" region. Figure lb shows the resultant spectra, where the 
background spectra created from eight (SL) and four (LL) slit-aperture data at both ends of 
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Table 1. Observation Log 



Instrument 


Date 


Observing mode 


Observation ID 


Details 


AKARI/IRC 


Nov 12, 2006 


IRC02 


1400041 


N3, N4, S7, Sll band imaging 


AKARI/IRC 


May 11, 2007 


IRC02 


1402057 


LI 5, L24 band imaging 


AKARI/FIS 


Nov 13, 2006 


FISOl 


1400044 


Scan speed: 8" s^^. Reset interval: 2.0 s 


Spitzer/IRS 


May 13, 2009 


IrsMap 


26091264 


see text 



each mapping area are subtracted. We observed a small mismatch between the continuum 
levels of the SL and LL spectra and hence applied a scale factor of 1.2 to the LL data to obtain 
their agreement at the overlapping wavelengths. 

The AKARI data were taken in part of the AKARI mission program "ISM in our Galaxy 
and Nearby Galaxies" (ISMGN; PI: H. K.). With the IRC, we obtained the N3 (reference 
wavelength of 3.2 ^m), N4 (4.1 yum), 57(7.0 ^m), Sll (11.0 /xm), L15 (15.0 ^m), and L24 (24.0 
lim) band images using a standard staring observation mode. The images were created by using 
the IRC imaging pipeline software version 20071017 (see IRC Data User Manual; Lorente et 
al. 2007 for details). The background levels were estimated by averaging values from multiple 
apertures placed around the galaxy, and subtracted from the images. With the FIS, we obtained 
N60 (centered at a wavelength of 65 ^m), WIDE-S (90 /xm), WIDE-L (140 /im), and N160 (160 
ixm) band images using a standard slow-scan observation mode, a 2-round-trip slow scan with 
a scan speed of 8" s~^. A region approximately 10' x 10' around the galaxy was covered with 
a scan map. The images were processed from the FIS Time Series Data using the AKARI 
official pipeline developed by the AKARI data reduction team (Verdugo et al. 2007) with 
custom procedures to remove cosmic-ray-induced detector artifacts. The background levels 
were estimated from data taken near the beginning and the end of the slow scan observation 
and subtracted from the images. 

3. Results 

3.1. Spitzer mid-IR spectral maps 

As seen in figure lb, we clearly detect the prominent PAH 11.3 fim feature and [Sill] 34.8 
/im emission line, together with several emission lines such as H2S(3) at 9.7 i^im, [Nell] at 12.8 
fj,m, [Nelll] at 15.6 /im, H2S(1) at 17.0 fim, [Sill] at 18.7 and 33.5 /im, the PAH 17 fim broad 
feature, and the mid-IR dust continuum emission at wavelengths longer than 27 fim. Figure 2 
shows the contour maps of the spectral components for the central 45" x 45" area, where the 
bin sizes of spectral images are set to be 1"85 for SL and 5''08 for LL. We applied smoothing 
with a boxcar kernel of 3 pixels in width (~ 5" for SL and ~ 15" for LL) for every spectral 
map. In table 2, we list the wavelength ranges used for extracting each spectral component. 
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Fig. 1. (a) illustrative views of IRS SL (left) and LL (right) slit positions for the spectral mapping obser- 
vations of NGC 4125, superposed on the 2MASS J-band image of NGC 4125. The image sizes are 1'.3 x 1'.3 
for SL and 3' x 3' for LL. The north is up and the east is to the left, (b) The background-subtracted IRS 
SL (blue) and LL (red) spectra of the central 15" x 15" region of NGC 4125. 
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Fig. 2. Spitcr/IRS spectral maps of the 5.5 — 6.5 /im continuum, 14 — 15 /im continuum, 30 — 35 /im 
continuum, H2S(3) 9.7 fim, PAH 11.3 fim, [Nell] 12.8 /im, PAH 17 fim, H2S(1) 17.0 /im, and [SiH] 34.8 
fim emissions. For each map, the contours are drawn on a hnear scale from 90 to 10 % of the peak surface 
brightness. The 5 % and 7 % level contours are added for the 5.5 — 6.5 /im continuum and PAH 11.3 /im 
map, while the 5 % level contour is added for the H2S(1) map. 

Note that the [Nell] and H2S(1) hne emissions are blended with the underlying PAH 12.7 fim 
and 17 fim features, respectively. For the PAH 6.3 and 7.7 /im features, we could not obtain 
reliable spectral maps because the slopes of the underlying continua are steep and the features 
are broad. 

As seen in figure 2, the distributions of the continuum emissions at wavelengths of 5.5- 
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Table 2. Wavelength ranges used for deriving the spectral maps in Fig. 2 



LyOmponent 


Signal 


Backg 


;round 




(/.tm) 




(/im) 


5.5 — 0.5 continuum 


5.5 — D.5 






14—15 fim continuum 


AAA 1 /I n 

14.4 — 14. y 






30 — So i-im continuum 


30.0 — 34.5 






TT O /'O \ A T 

H2b(oj 9.7 /im 


A AO 

9.0 — 9.8 


An AC An 

9.0 — 9.5, 9.9 


1 A A 

— 10.9 


rAU ii.o /ini 


11.0 — 11.0 


An iin 110 
y.9 — 11.0, 11.(3 


— lz.3 


[NcII] 12.8 ^lm. 


12.8-12.9 


13.0 


- 13.2 


PAH 17 


16.5-16.9, 17.3-18.6 


14.4-15.5, 19.5 


-20.7 


H2S(1) 17.0 fim 


16.9-17.3 


14.4-15.4, 15.8-16.3, 19.4 


-20.8 


[Sill] 34.8 /im 


34.5-35.5 


33.8-34.5, 35.5 


-36.2 



6.5 /im and 14-15 //m show an elliptical shape, which is consistent with the stellar distribution 
in the 2MASS image (figure la), while that at 30-35 //m is of a more circular shape. These 
distributions support that the continuum emission at wavelengths shorter than 20 /im is dom- 
inated by photospheric emission from old stars as also suggested by the spectral slope, and 
that the mid-IR dust continuum emission at wavelengths longer than 27 /im is distributed dif- 
ferently from the stellar emission. In fact it is generally known that near- to mid-IR emission 
from an elliptical galaxy is dominated by photospheric emission from low-mass giant stars with 
temperatures of 3000 K to 5000 K (Athey et al. 2002; Xilouris et al. 2004; Temi et al. 2008). 

The distributions of the PAH 11.3 /im and 17 /im as well as those of H2S(3) and S(l) 
are relatively compact as compared to the stellar distributions at 5.5-6.5 /im and 14-15 /xm. 
The PAHs and H2 are extended in a similar direction with the position angle (P.A.) of ~90°, 
where P.A. is measured from the north through the east. The atomic gas line emissions of 
[Nell] (at relatively high brightness levels) and [Sill] are extended similarly in the direction of 
P.A. of ~70°, but differently from the PAHs and H2. More interestingly, they are extended 
in different directions from the major axis P.A. 82.5° of the stellar distribution (Jarrett et al. 
2003). Quantitative comparisons of spatial distributions among all the relevant images are 
summarized later in the discussion (section 4.3). 

In order to clearly visualize the difference in size of the distribution between each spectral 
component, spectra are compared in figure 3, which are created by integrating the areas of 
5" X 5", 10" X 10", 15" X 15", and 25" x 25". As seen in figure 3a, the intensities of some spectral 
features with respect to the underlying continuum intensities change from the inner to the outer 
area. Figure 3b shows the fiux ratios of the innermost (5" x 5") to the outermost (25" x 25") 
spectrum. Since the smallest area is comparable to the beam size of the IRS/LL module, we 
cannot discuss the change of the continuum shape. But the change of the PAH features and the 
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line emissions relative to the underlying continua (i.e. their equivalent widths) can be discussed 
from figure 3b. It is clear that only the PAH 11.3 fim and H2S(3) emissions, and possibly the 
H2S(1) emission change significantly, while all the atomic gas line emissions show no significant 
changes, implying that the PAHs and H2 are the most centrally concentrated. In particular, 
the difference between the PAH 11.3 /im and [Nell] 12.8 /im emission is conspicuous; there is 
no change in the [Nell] line equivalent width from the inner to the outer area, demonstrating 
that the distribution of the [Nell] emission is extended more than the PAH 11.3 fim emission. 
The structures around 14 /im and 19 fim are probably not real because they are located near 
the edges of the 1st order of the SL module and the 2nd order of the LL module, respectively. 

3.2. AKARI near- to far-IR images 

The near- to far-IR 9-band images of NGC 4125 are shown in figure 4. We do not use 
the N160 data in the following discussion because we could not obtain significant detection. 
For every map in figure 4, the image size is set to be the same (~ 3' x 3'). The bin sizes 
of the images are IM6 for N3 and A^^, 2:'34 for S7 and Sll, 2'.% for L15 and 7'.'5 for 
N60 and WIDE-S, and 12"5 for WIDE-L. As seen in figure 4, the spatial distribution changes 
considerably with the photometric band. The N3, N4 and S7 images show smooth distributions 
tracing the stellar component of the galaxy. However, the Sll and Li5 bands start to deviate 
from the smooth stellar distribution. The L24 band image shows a relatively faint compact 
emission. These changes of the images in the mid-IR are quite reasonable judging from the 
above Spitzer spectra. As for the 11 /im-bright spot located at ~ V to the southwest from the 
center, there is no counterpart found in the SIMBAD database. This can also be recognized in 
the 15 /im and 24 fim bands and hence is probably a real mid-IR source. However the 11 fim 
emission bridging the gap between the galaxy and the source is probably caused by a low-level 
ghost signal of the galaxy emission due to the internal refiection within the beam splitter, which 
are known to appear at about 1' from an original signal. We ignore this southwest emission 
structure below. 

The N60, WIDE-S, and WIDE-L images, again, show apparent deviations from the 
smooth stellar distribution as seen in the N3 and N4 bands. The WIDE-S and WIDE-L images 
exhibit an apparent elongation to a similar direction to the [Nell] and [Sill] line maps in figure 2 
(i.e. P. A. ~ 70°), although their spatial scales are quite different. As explained in Kaneda et al. 
(2008b), the WIDE-S band image is best among the three and even advantageous over MIPS 
far-IR images to discuss the spatial distribution of faint far-IR emission in elliptical galaxies. 
In figure 5, we show the encircled energy profile of the far-IR distribution of NGC 4125 in 
the WIDE-S band, comparing to that of the point spread function (PSF) in the same band 
(Shirahata et al. 2009). From the figure, we conclude that the distribution of the far-IR dust 
emission is significantly extended even with the large PSF of the far-IR band. For comparison, 
we also show an encircled energy profile for the PAH 11.3 fim emission in figure 2. It is obvious 
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Fig. 3. (a) Spitzer/IRS spectra created by integrating data witliin tlie central 5" x 5", 10" x 10", 15" x 15", 
and 25" x 25" areas in the ascending order, (b) The ratio of the spectrum of 5" x 5" to that of 25" x 25". 
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Fig. 4. AKARI/FIS and IRC images of NGC 4125. The contours arc drawn at logarithmically-spaced 10 
levels from 80 % to 2 % of the peak surface brightness. The contour levels below 3 %, 5 %, and 10 % are 
not shown for the L24^ and the FIS images, respectively. 
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Fig. 5. Encircled energy profile of the far-IR distribution of NGC 4125 in the WIDE-S band in open 
triangles, which is compared to that of the standard PSF in the same band in solid curve (Shirahata et 
al. 2009). The two profiles are normalized to unity at a radius of 180". The open squares indicate the 
encircled energy profile of the distribution of the PAH 11.3 /im emission. The error bars attached to the 
open triangles and squares are estimated from sky fluctuations in their image data, while those to the 
solid curve are the standard deviations of the observed PSFs with different colors (Shirahata et al. 2009). 

that the distribution of the PAH emission is by far more compact than the PSF in the far-IR, 
and that there is a large difference in size of the distribution between the PAHs and far-IR 
dust. 

3.3. Spectral energy distributions and spectral fitting 

We derive the flux densities of NGC 4125 by integrating the surface brightness within 
the circular aperture radii of 14f'l (~ 1.5 kpc) and 1^5 (~ 9.7 kpc) around the galactic center 
in each AKARI band image. The former aperture has the same area as the largest one for 
the Spitzer/IRS specta (figure 3). The observed flux densities are listed in table 3, where the 
errors are estimated from background sky fluctuations. For the IRC, no aperture corrections 
are performed, since the aperture radii are sufficiently large as compared to the PSF sizes. 
For the FIS, aperture corrections are performed by using the correction tables in Shirahata et 
al. (2009), while color corrections are not performed because they are found to be negligible 
considering the systematic errors of the fluxes described below. 

The spectral energy distributions (SEDs) constructed from the flux densities in table 3 
are presented in flgure 6. For the SED of the central region (r < 14f'l), the above Spitzer/IRS 
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spectrum of the same area is shown together, while for the SED of the total region (r < 1'5), 
the Spitzer/IRAC and MIPS photometric data points are added (Dale et al. 2007; Temi et al. 
2007). For the FIS data points, we consider systematic errors associated with the detectors 
and absolute uncertainties, which are estimated to be 20 % for N60 and WIDE-S, and 30 % for 
WIDE-L (Kawada et al. 2007). The AKARI flux densities in the near- to far-IR show overall 
agreements with the measurements by Spitzer including both photometric data with the IRAC 
and MIPS and spectral data with the IRS. 

In order to fit the IRS spectra, we have used the PAHFIT tool developed by the Spitzer 
SINGS legacy team (Smith et al. 2007b). In table 4, the fitting results of the spectra for the 
15" X 15" and 25" x 25" areas are summarized for the intensities of the gas lines, PAH features 
and continuum emissions that have signal-to-noise ratios higher than 5 for both spectra. To 
further fit the far-IR photometric data points together with the mid-IR spectra in figure 6, we 
have changed the temperature of the coolest dust component in the PAHFIT from 35 K to 25 
K (figure 6a) or 23 K (figure 6b); all the dust components have the emissivitiy power-law index 
of /3 = 2. We have simply extrapolated the best-fit curve in the mid-IR to the near-IR regime. 
For the total SED in figure 6b, the mid-IR spectral shape is kept the same as the above best-fit 
curve and only the normalizations are scaled to match the photometry. 

The model reproduces both SEDs well. The Sll data point exceeds the level interpolated 
by using the other photometric data points, reflecting that the strong PAH 11.3 /xm feature 
and its underlying plateau emission contribute much to the SI 1 band intensity. The S7 band 
intensity is dominated by the stellar continuum emission like the N3 and the band intensity, 
because usually strong PAH 7.7 /im band emission is apparently weak; the PAHFIT gives the 
PAH 7.7 //m complex intensity of (5.9 ±1.2) x 10"^'' W m~^ for the 15" x 15" area and no 
significant detection (< 1.2 x 10^^^ W m~^ as a 5-sigma upper limit) for the 25" x 25" area. 
It should be noted, however, that the PAH 7.7 /im band strength is underestimated unless 
a stellar silicate feature around 10 yum is taken into account, and that the difference can be 
such a large factor as 2-3 for the above feature strength (Kaneda et al. 2008a). Nevertheless, 
even considering the effect of the stellar component, the PAH7.7/11.3 ratio is 0.6-0.9, which 
is considerably lower than 3.6, the median, and 1.5-4.8, the 10 %-90 % range of variation 
over the full sample of the normal galaxies in the Spitzer/SINGS legacy program (Smith et al. 
2007b). The PAH 17/11.3 ratio is 0.5 (table 4), very similar to the ratio of ~0.5 averaged over 
all the SINGS sample, suggesting that the size distribution of the PAHs is normal (Draine & 
Li 2007). Thus we attribute this unusual weakness of the PAH 7.7 /im complex intensity to the 
dominance of neutral PAHs over ionized ones due to very soft interstellar radiation field typical 
of elliptical galaxies (Kaneda et al. 2005, 2008a). However we detect the PAH 6.2 //m feature 
that is apparently broad and strong despite the weak PAH 7.7 /im feature. Consequently, the 
PAH6.2/11.3 ratio is 0.8-1.1 (table 3), similar to 1.1 and 0.52-1.5 for the above values of the 
SINGS sample. By referring to the diagnostic diagram of PAH 6.2 /im/7.7 /im versus 11.3 
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Table 3. Observed Flux Densities of NGC 4125 



Band name 


Arcf* 


Center (r < 14"1) 


Total (r < 1.'5) 


IRC N3 


3.2 fini : 


197±5 mJy 


574±15 mJy 


IRC N4 


4.1 /ini : 


117±4 mJy 


339±11 mJy 


IRC 57 


7.0 /ini : 


61±2 mJy 


178±4 mJy 


IRC Sll 


11.0 /ini : 


49±1 mJy 


145±3 mJy 


IRC L15 


15.0 fiia : 


32±1 mJy 


105±3 mJy 


IRC L2J^ 


24.0 : 


26±2 mJy 


80±5 mJy 


FIS N60 


65 fiia : 


361±47 mJy 


575±34 mJy 


FIS WIDE-S 


90 /ini : 


525±58 mJy 


953±38 mJy 


FIS WIDE-L 


140 /im : 


573±139 mJy 


1250±83 mJy 



* Reference wavelengths (Onaka et al. 2007; Kawada et al. 2007) 

/im/7.7 fim ratios in Sales et al. (2010), our result reveals that the 6.2 ^m/7.7 fim ratio is 
remarkably large. This might suggest the relative importance of substitution of C by impurity 
atoms (Hudgins et al. 2005) or clustering of PAHs for the PAH 6.2 /xm feature (Rapacioli et 
al. 2005; Tielens 2008), although we do not find any other hints from the other PAH features. 

3.4- Physical parameters derived from spectral fitting 

First we obtain a robust estimate on the far-IR (50-170 /im) luminosity, Lakari? since 
the FIS N60, WIDE-S, and WIDE-L bands continuously cover the wavelength range from 50 
/xm to 170 /im (Hirashita et al. 2007). The results are Lakari = 2.8 x lO^L© for r < 14''1 and 
i^AKARi = 5.1 X IO^Lq for r < 1^5. 

We calculate the dust mass, Mdust, by using the equation (e.g. Hildebrand 1983): 

Mlust = „ /^x, (1) 

where D, F^, n^, and By{T) are the galaxy distance, the observed flux density, a dust mass 
absorption coefficient, and the value of the Planck function at the frequency of v and the dust 
temperature of T, respectively. We use the 90 /xm flux density and the dust temperature, both 
given by each dust component in figure 6. We adopt a dust mass absorption coefficient of 36 
cm^g"^ at 90 /tm given by Hildebrand (1983) and adjusted for (3 = 2. The dust mass thus 
derived is 3.4 x 10^ Mq for r < 14"1 and 1.0 x 10^ Mq for r <1%. Considering the difference in 
adopted galaxy distance, the latter dust mass is ~ 3 times larger than that derived with IRAS 
by Goudfrooij & de Jong (1995), while it is ~ 2.5 times smaller than the mass derived with 
ISO by Temi et al. (2004) who included a cooler (20 K) dust component. Even colder dust is, 
however, unlikely to exist predominantly because the IRAM observation of NGC 4125 at 1250 
/xm gave 3.5±1.9 mJy, only a tentative detection (Wiklind & Henkel 1995). Thus the absolute 
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Fig. 6. (a) AKARI photometric data points (boxes with the error bars) for a circular integration aperture 
of r ~ 14"! together with the Spitzcr spectral data points (with the dots) for the same integration area. 
The error bars in the wavelength indicate the effective band widths. The green curve represents the curve 
fitted by the PAHFIT program with a newly-added cooler dust component. The contribution by dust 
gray-body components with different temperatures is shown by the red curves, while the stellar emission 
component is shown by the magenta curve, (b) AKARI photometric data points for a circular integration 
aperture of r = 1.'5 fitted by the above model where the mid-IR spectral shape is kept the same as in (a) 
but only the normalization is scaled to match the photometry. The crosses with error bars only in the flux 
correspond to the Spitzcr photometric data points (Dale et al. 2007; Temi et al. 2007). 
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values of the dust masses are quite uncertain, but the relative change from the inner to the 
outer region suggests that a significant fraction of dust is distributed over the galaxy. 

From the pair of the ortho-H2 S(l) and S(3) line intensities in table 4, we calculate the 
temperature and column density of the molecular hydrogen gas by using an excitation diagram 
with the transition probabilities of the H2 pure rotational lines from Black & Dalgarno (1976). 
The results are T = 300 K and Nu, = 4.1 x 10^^ cm'^ for the 15" x 15" area and T = 310 K and 
A^Ha = 1-7 X 10^* cm"^ for the 25" x 25" area. Thus the warm H2 masses are 1.7 x 10^ Mq for 
the former and 2.0 x 10^ Mq for the latter, which are comparable to the dust mass. The cold 
H2 mass is estimated to be 7 x 10'^ Mq from tentative CO detection by Wiklind et al. (1995), 
which gives the gas-to-dust mass ratio of 70-200, similar to the accepted value of 100-200 for 
our Galaxy (Sodroski et al. 1997). 

On the other hand, the line pair of the [Sill] 33 /xm and 19 /xm lines is a good measure 
of an electron density rig of ionized gas. From table 4, their intensity ratios correspond to 
Ue ~ 350 for the 15" x 15" area and < 30 for the 25" x 25" area (Rubin et al. 1994). The 
highly-ionized gas responsible for the lines such as [Sill], [Nell], and [Nelll] is likely to be 
associated with the LINER nucleus (Smith et al. 2007b). Alternatively, the highly-ionized gas 
emission might be powered by interaction with the hot plasma through heat conduction (Sparks 
et al. 1989). Because the [Nell] line emission is extended in similar directions to the X-ray 
plasma (see below), we prefer the latter scenario. As clearly demonstrated in figure 3b, the 
[Nell] line emission is significantly more extended than the PAH 11.3 fim emission. 

As for the relative intensities of the dust, molecular, and atomic gas emissions, the ratio 
of the [Sill] line intensity to the total AKARI far-IR flux is 4.4 x 10~^ for the aperture of the 
same area (25" x 25" or r < 14"1), which is considerably higher than 2 x 10""^, the average of 
the SINGS sample galaxies (Roussel et al. 2007). The ratio of the H2 lines (S(l)-S(3)) to the 
[Sill] line intensity for the central 15" x 15" area is also as high as ~1.5, which is well outside 
the range of 0.15-0.5 for the nuclei of the SINGS sample galaxies (Roussel et al. 2007). Hence 
the H2 line emissions are extremely bright and the [Sill] line is also strong relative to the far-IR 
emission, as compared to the canonical relation for the normal galaxies, suggesting that they do 
not originate simply from photo-dissociated regions. There might be a dominant contribution 
from slow shock heating by dynamical perturbation of gas clouds in the center for the former 
(Roussel et al. 2007) and Si enrichment through sputtering destruction of dust by hot plasma 
for the latter (Kaneda et al. 2008a). The strong [Sill] line may also be due to the cooling of 
X-ray dominated regions around the LINER nucleus (Dale et al. 2006). 
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Table 4. Intensities of gas lines, PAH features, and continuum emissions 



Are& 


[Nell] 


[Nelll] 


[SIII] 


fSIIIl 


[Sill] 




H9Sf3~) 




12 82 um 


15.55 /iin 


18.71 fim 


33 45 ;;m 


34 77 iim 


17 04 /vm 


Q fifi //TD 


15" X 15" 


2 36±0 in* 


2.04±0.08 


0.88±0.09 


1.11±0.12 


4 55±0 18 


3 24±n 1 5 


2 28±0 09 


25" X 25" 


3.08±0.17 


2.52±0.21 


1.23±0.15 


2.35±0.24 


7.93±0.26 


4.00±0.26 


3.35±0.15 


Area 




PAH 






Continuum 






6.2 (Um 


11.3 /im 


12.7 


17 /im 


6 /im 


15 /im 


34 /im 


15" X 15" 


16.3±1.2 


19.81±0.39 


4.61±0.85 


10.39±0.66 


42.74±0.33 


15.49±0.45 


31.4±1.4 


25" X 25" 


26.5±2.3 


24.25±0.63 


5.77±0.50 


12.7±1.3 


68.27±0.53 


27.3±1.1 


57.3±3.3 



* The values are given in units of mjy for the continuum emissions and 10 W m ^ for the feature and the line emissions. 

4. Discussion 

4.I. PAHs in the central dense gas region 

We discuss the distribution of the PAH emission. In figure 7a, the 5.5 — 6.5 /im continuum 
emission shows a smooth stellar distribution, whereas the PAH 11.3 /im emission exhibits an 
apparently compact and different distribution from the stars. Since PAHs can be excited even 
by visible and near-IR soft photons from old stars (Li & Draine 2002), the compact distribution 
of the PAH 11.3 /im emission strongly suggests the real absence of neutral PAHs in the outer 
regions of the galaxy. It is possible that PAHs are increasingly becoming charged and the 11.3 
/im emission becomes weaker from the center to the outward (Draine & Li 2001; Li & Draine 
2001). However the PAH 7.7 /xm emission, which traces the distribution of charged PAHs, does 
not tend to become stronger in the outer regions, as seen in Fig. 3a. A more likely interpretation 
is that PAHs have been destroyed in the X-ray hot plasma outside the central dense gas region, 
because PAHs are expected to be destroyed in so short a timescale as 10^ yr for the plasma 
parameters typical of X-ray elliptical galaxies (Micelotta et al. 2010; see below). Deep in dense 
molecular gas, PAHs are shielded and may be efficiently hydrogenated, which can explain the 
observed strong 11.3 /im emission. 

Figures 7b and 7c show comparisons of the PAH emission with the V — I color in- 
dex image and the continuum-subtracted Ha image, respectively, which are taken from the 
NASA/IPAC Extragalactic Database. The former image reveals the presence of dense gas re- 
gions near the galactic center, and the PAH emission follows its distribution. Verdoes Kleijn 
& de Zeeuw (2005) observed a central dust lane of size lf'3 with the difference in P. A. of +13° 
from the galaxy major axis. Although the spatial scale is much larger than this, the PAH 
emission region has a spatial correspondence with the optical dust lane similarly to the case 
of NGC 4589, where the PAH emission comes predominantly from the optical dust lane of the 
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elliptical galaxy (Kaneda et al. 2010). It is notable that the distribution of the Ha emission 
also resembles those of the dust lane and PAHs. It is likely that the Ha emission is not simply 
related to star formation but rather cooling gas from the hot plasma or heating gas from the 
dense cold gas through thermal conductance at their contacting interfaces. 

As seen in the spectra in figure 6, the Sll band is expected to contain strong PAH 
11.3 /im emission and its underlying plateau feature probably due to PAH clusters, in addition 
to the stellar continuum emission. The contribution of the [Nell] line is negligible from table 
4. Since the S7 band is dominated by the stellar continuum emission, we can estimate the 
spatial distributions of the PAHs and PAH clusters emitting the features around 11 /im by 
normalizing the 5*7 image properly and subtracting it from the Sll image. We determined the 
normalization factor on the basis of the stellar continuum predicted by the spectral fitting in 
figure 6a. In the S'ii- 5*7 differential image thus derived (figure 7d), the emission extends along 
the east-west direction (i.e. P.A.= 90°) rather differently from the major axis P.A. 82.5° of 
the stellar distribution, and connects to the inner distribution of the PAH 11.3 /im emission in 
figure 7a. The broader distribution than the PAH 11.3 //m emission is probably attributed to 
the PAH clusters responsible for the 11 /im plateau emission. 

4-2. Dust with diffuse hot plasma and atomic gas 

We compare the distribution of dust with those of the X-ray plasma and the atomic 
gas. Figure 8a shows the X-ray (0.3-7 keV) image retrieved from the Chandra data archive. 
In order to obtain the spatial distribution of diffuse hot plasma alone. X-ray point sources that 
are thought to be low-mass X-ray binaries were removed from the image by using the Chandra 
Interactive Analysis of Observations (CIAO) software. The X-ray image was then smoothed and 
the contours were superposed on the original image. As seen in figure 8a, NGC 4125 certainly 
possesses diffuse hot plasma, but its distribution is rather different from the stellar distribution. 
Figure 8b shows a comparison between the X-ray plasma and the stellar continuum emission in 
the central 45" x 45" area, revealing that the X-ray distribution has a smaller P.A (~70°) than 
the stellar major axis P.A. (~ 82.°5) of the galaxy. Figures 8c and 8d show the X-ray images 
in a soft (0.3 — 0.8 keV) and a hard (0.8-7 keV) band, respectively, which clearly exhibit that 
cooler plasma is distributed more widely to the direction of P.A.~70°. Trinchieri et al. (2000) 
already pointed out that the ASCA and BeppoSAX 0.2-10 keV spectrum of NGC 4125 cannot 
be reproduced by a single thermal model, calling for the presence of a second soft (kT~0.3-0.7 
keV) spectral component. 

Figures 9a and 9b show the spatial distribution of the far-IR dust compared to those 
of the atomic gas ([Sill]) and X-ray hot plasma, respectively; they bear resemblances to each 
other, which are extended to similar directions of P.A.~70°. The spatial correlation in figure 
9a may suggest relatively abundant gas-phase Si through plasma destruction of dust where Si 
is depleted. As shown in figure 9b, the X-ray hot plamsa is obviously of a diffuse nature, and 
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Fig. 7. (a) Spitzcr/IRS contour map of the PAH 11.3 /xm feature (red lines) overlaid on that of the 
5.5 — 6.5 /xm continuum emission (black lines). Both are the same as in figure 2. (b) The V — I color index 
image superposed on the same contour maps as in the panel (a), (c) The continuum-subtracted Ha image 
superposed on the same PAH 11.3 /im map as in the panel (a), (d) AKARI/IRC Sll image shown after 
subtraction of the S7 image with its intensity adjusted on the basis of the SED in figure 6a. The contours 
are drawn at logarithmically-spaced 10 levels from 80 % to 2 % of the peak surface brightness. Note that 
spatial scales are different between the panels (a)~(c) and the panel (d) 
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the dust emission is distributed similarly; such a relationship is rather unexpected because the 
lifetime of dust with sizes of 0.1 /xm exposed to hot plasma with gas densities of 10~^ ~ 10"^ 
cm~^ is quite short (10^ ~ 10^ yr; Draine & Salpeter 1979; Tielens et al. 1994; Tsai & Mathews 
1995). NGC 4125 has electron densities and temperatures of (2 — 10) x 10^^ cm"^ and 0.36-0.44 
keV at 3-10 kpc radii (Fukazawa et al. 2006). The supply of the dust by old stars as dominant 
sources is unlikely considering the clear difference in distribution between the dust and stars. 
If the observed dust mostly originates from external sources, their lifetime is obviously shorter 
than the crossing time of a galaxy through a merging galaxy (~ 10* yr), where we assume the 
galaxy size of 10 kpc diameter with an infall velocity of 100 km s~^. A possible explanation 
is the evaporation flow scenario (Sparks et al. 1989; de Jong et al. 1990), where evaporating 
clouds were brought in by a gas- rich galaxy that merged in the past. The evaporation of the 
clouds is due to thermal conduction by plasma eletrons resulting in cooling the X-ray plasma. 
For NGC 4696, which has X-ray plasma parameters similar to the above, de Jong et al. (1990) 
calculated that clouds with column densities of iVn ~ 7 x 10^^ cm~^ or Ay — 3 can survive in 
10® yr. The alignment of the distributions of the dust to the softer plasma emission supports 
the evaporation cloud scenario. In the present case, however, clouds must be sufficiently small, 
and thus dense, because they are not observed in the extinction map (figure 7b). Even if clouds 
escape detection in the optical absorption, they would have to be observed in the PAH emission. 
However, the observed PAH emission does not have such extended components, either. In fact, 
the ratio of the integrated PAH (6.2 — 17 ^m) intensity to the total far-IR fiux is 3.5 x 10"^ 
for the area of 25" x 25" (or r < 14"1), which is considerably lower than 0.10 ±0.04 for the full 
sample of SINGS galaxies (Smith et al. 2007b). Therefore it is difficult to explain the observed 
dust emission by the presence of such dense clumpy clouds in the outer region of the galaxy. 

4.3. Constraints on origins of dust and PAHs from their distributions 

Figure 10 summarizes the size and the axis direction of the spatial distribution of each 
spectral and photometric component in terms of the outer/inner flux ratio and the P.A. The size 
of the distribution is normalized to that for the 4 /im band intensity, i.e. the stellar emission. 
It is clear from the flgure that the emissions by the PAHs as well as the molecular and atomic 
gases excluding the [Sill] emission are more centrally concentrated than the stellar emission, 
and on a larger scale, the emissions by the far-IR dust and the X-ray hot plasma are not more 
widely distributed than the stellar emission (probably less for the far-IR emission). All of these 
suggest that the PAHs and dust are not of an external origin; otherwise their distributions 
would be found to be extended more widely than the stellar distribution. In flgure 10, we flnd 
that there are at least three types of axis directions. The continuum and the photometric band 
emissions dominated by the stellar emission possess the P.A. of 82.°5, while the atomic gas, 
dust, and hot plasma show ~ 70°, a significanly smaller P.A. The PAHs and molecular gas 
show ~ 90°, a significanly larger P.A. than the stellar emission. 
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Fig. 8. Chandra X-ray smoothed image of NGC 4125. (a) The 0.3 — 7 keV X-ray contour map after 
removing point sources and smoothing an image, overlaid on the original X-ray image, (b) The same 
contour map as in (a) on that of the 5.5 — 6.5 yum continuum emission, (c) The 0.3 — 0.8 keV and 0.8 — 7 
keV X-ray contour maps. The X-ray contours are drawn at logarithmically-spaced 10 levels from 80 % to 
3 % for the panels (a) and (b) and from 80 % to 10 % for the panels (c) and (d). 



20 




Fig. 9. (a) contour map of the [Sill] 34.8 fim emission, the same as in figure 2, superposed on the WIDE-S 
90 fim image, the same as in figure 4, but displayed in the color levels, (b) The same X-ray contour map 
as in Fig. 8(a) superposed on the 90 ^m image. Note that the image data of [Sill] are limited only to the 
central 45" x 45" area. 

Near the center of the galaxy, the PAHs and dust are hkely to have been protected 
from the interaction with the X-ray plasma inside the dense molecular gas; otherwise the PAHs 
would have been destroyed in a very short time scale of ~100 yr (Micelotta et al. 2010). In 
other words, the hot plasma does not occupy the whole interstellar space near the center of 
the galaxy. In the outer region of the galaxy, however, from the similarity in their spatial 
distributions, the dust and hot plasma must be mixed to some extent, unless the volume filling 
factor of the hot plasma is small. It is very unlikely that the hot plasma has non-diffuse patchy 
distributions. In fact, our observation suggests that the dust is interacting with the hot plasma, 
producing cooler plasmas. As described above, dust grains with a 0.1 fim size can survive for 
a 10^ — 10'' yr timescale, which is much longer than the lifetime of the PAHs but much shorter 
than the look-back time of (4 — 6) x 10^ yr for the merger event (Schweizer & Seitzer 1992). 
Therefore the dust is likely to have been supplied into the outer region of the galaxy quite 
recently apart from the merger event. 

Our observation shows that PAHs are present only in the innermost region while the 
dust emission is obviously extended than the PAH emission. If we regard the observed PAHs as 
smallest form of dust grains, the observational result suggests the size- dependent distribution 
of dust, i.e. larger grains extended to outer areas, which favors an internal outflow origin for 
the dust-lane and extended materials. Temi et al. (2007) suggested that the buoyant outflow 
of dust from a central reservoir in the core of a galaxy may play an important role in supplying 
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Fig. 10. (Top) Ratios of fluxes in outer to inner area for various spectral components and photometric 
band intensities, which are normalized to the corresponding ratios in the N4 band. The triangles represent 
the ratios in the outer 25" x 25" to the inner 15" x 15" area, while the boxes are those in the r = 
area to the inner r = 14"1 area. (Bottom) The position angles of the spatial distributions for the various 
components obtained by an ellipse fitting. The triangles and the boxes represent the results of fitting to 
15 % and 20 % brightness levels of the peaks, respectively. Only for the X-ray image arc they 20 % and 
30 % brightness levels. The errors are obtained by assuming the positional uncertainty of 0.5 pixels per 
bin. The dashed line corresponds to the stellar major axis of the galaxy (P.A.= 82.°5; Jarrett et al. 2003). 



dust into the interstellar space of elliptical galaxies with a rather high buoyant velocity up to 
400 km s~^ possibly by receiving momentum from a jet emerging from a central nucleus. From 
this viewpoint, the extended dust and cooling gas emission structure might be related to the 
outflows from the LINER dusty nucleus of NGC 4125. 

5. Conclusions 

We have performed mid-IR spectral mapping with Spitzer and near- to far-IR imaging 
observations with AKARI for NGC 4125, a nearby elliptical galaxy possessing diffuse interstellar 
X-ray hot plasma. We obtain the spatial distributions of dust and PAHs as well as molecular and 
atomic gases in NGC 4125. We also derive the distribution of diffuse hot plasma by removing 
point-like sources from Chandra X-ray images, and compare them to understand how they are 
separated or coexist in the interstellar space. We find that the spatial distributions of the 
PAHs and dust are different from each other. The PAH emission predominantly comes from a 
compact dense molecular gas region near the center of the galaxy, and hence the observed PAHs 
are likely to have been protected from the interaction with the X-ray plasma inside the dense 
molecular gas. There are no diffuse components observed in the PAH emission. In contrast, the 
dust and atomic gas emissions have more extended structures similar to the distribution of the 
X-ray plasma. The similarity suggests their association with the X-ray plasma. The PAH and 
the dust emission do not follow the stellar surface brightness, suggesting unimportance of stellar 
mass loss as their sources. NGC 4125 is thought to be a (6 — 8) x 10^ yr old remnant of a major 
disk-disk merger. Even if dust and PAHs were brought in by the merger, they cannot survive 
until now considering their short destruction timescales in hot plasma. The PAHs and dust 
are likely to have been produced secondarily after the merger, having avoided their interaction 
with the hot plasma for a long time near the center of the galaxy. Recent transient events such 
as outflows from the central dusty nucleus may have spread the dust over the outer region of 
the galaxy. 

We would express many thanks to an anonymous referee for giving us useful comments. 
This work is based on observations made with AKARI and Spitzer. AKARI is a JAXA project 
with the participation of ESA. Spitzer is operated by the Jet Propulsion Laboratory, California 
Institute of Technology under NASA contract 1407. We are grateful to the SINGS team for 
providing the CUBISM software and the PAHFIT tool. This work has also made use of data 
obtained from the Chandra Data Archive. This research is supported by the Grants-in-Aid for 
the scientific research No. 22340043, and the Nagoya University Global COE Program, "Quest 
for Fundamental Principles in the Universe", from the Ministry of Education, Culture, Sports, 
Science and Technology of Japan. 



23 



References 



Athey, A., Bregman, J. N., Bregman, J. D., Temi, P., & Sauvage, M. 2002, ApJ, 571, 272 
Bertola, F., Bettoni, D., Rusconi, L., k Sedmak, G. 1984, A J, 89, 356 
Black, J. H., & Dalgarno, A. 1976, ApJ, 203, 132 
Dale, D. A., et al. 2006, ApJ, 646, 161 
Dale, D. A., et al. 2007, ApJ, 655, 863 

de Jong, T., N0rgaard-Nielsen, H. U., J0rgensen, H. E., & Hansen, L. 1990, A&A, 232, 317 

Draine, B. T. & Li, A. 2001, ApJ, 551, 807 

Draine, B. T. & Li, A. 2007, ApJ, 657, 810 

Draine, B. T. & Salpeter, E. 1979, ApJ, 231, 77 

Fabbiano, G. & Schweizer, F. 1995, ApJ, 447, 572 

Filho, M. E., et al. 2002, ApJS, 142, 223 

Fukazawa, Y., Botoya-Nonesa, J. G., Pu, J., Ohto, A., &; Kawano, N. 2006, ApJ, 636, 698 

Goudfrooij, P. k de Jong, T. 1995, A&A, 298, 784 

Heckman, T. M., Balick, B., k Crane, P. C. 1980, A&AS, 40, 295 

Hildebrand, R. H. 1983, QJRAS, 24, 267 

Hirashita, H., et al. 2007, PAS J, 60, S477 

Houck, J. R., et al. 2004, ApJS, 154, 18 

Hudgins, D. M., Bauschlicher, C. W., Jr., k Allamandola, L. J. 2005, ApJ, 632, 316 
Humphrey, P. J., et al. 2006, ApJ, 646, 899 

Jarrett, T. H., Chester, T., Cutri, R., Schneider, S. E., & Huchra, J. P. 2003, AJ, 125, 525 

Jensen, J. B., Tonry, J. L., Barris, B. J., Thompson, R. I., Liu, M. C, Rieke, M. J., Ajhar, E. A., k 

Blakeslee, J. P. 2003, ApJ, 583, 712 
Kaneda, H., Onaka, T., k Sakon, L 2005, ApJ, 632, L83 

Kaneda, H., Onaka, T., Sakon, I., Kitayama, T., Okada, Y., k Suzuki, T. 2008a, ApJ, 684, 270 

Kaneda, H., Suzuki, T., Onaka, T., Okada, Y., k Sakon, I. 2008b, PASJ, 60, S467 

Kaneda, H., Onaka, T., Sakon, I., Kitayama, T., Okada, Y., Suzuki, T., Ishihara, D., k Yamagishi, 

M. 2010, ApJ, 716, L161 
Kawada, M., et al. 2007, PASJ, 59, S389 
Kennicutt, R. C, et al. 2003, PASP, 115, 928 

Knapp, G. R., Guhathakurta, P., Kim, D.-W., k Jura, M. 1989, ApJS, 70, 329 
Li, A. k Draine, B. T. 2001, ApJ, 554, 778 
Li, A. k Draine, B. T. 2002, ApJ, 572, 232 

Lorente, R., Onaka, T., Ita, Y., Ohyama, Y., k Pearson, C. P. 2007, AKARI IRC Data User Manual 
Version 1.2 

Micelotta, E. R., Jones, A. P., & Tielens, A. G. G. M. 2010, A&A, 510, A37 
Murakami, H., et al. 2007, PASJ, 59, S369 
Onaka, T., et al. 2007, PASJ, 59, S401 

O'Sullivan, E., Forbes, D. A., k Ponman, T. J. 2001, MNRAS, 328, 461 

Pu, S. B., Saglia, R. P., Fabricius, M. H., Thomas, J., Bender, R., k Han, Z. 2010, A&A, 516, A4 
Rapacioh, M., Joblin, C, k Boissel, P. 2005, A&A, 429, 193 

24 



Roussel, H., et al. 2007, ApJ, 669, 959 

Rubin, R. H., Simpson, J. P., Lord, S. D., Colgan, S. W., Erickson, E. F., & Haas, M. R. 1994, ApJ, 
420, 772 

Sales, D. A., Pastoriza, M. G., & Riffel, R. 2010, ApJ, 725, 605 
Schweizer, F. & Seitzer P. 1992, AJ, 104, 1039 
Shirahata, M., et al. 2009, PASJ, 61, 737 
Smith, J. D. T., et al. 2007a, PASP, 119, 1133 
Smith, J. D. T., et al. 2007b, ApJ, 656, 770 

Sodroski, T. J., Odegard, N., Arendt, R. G., Dwek, E., Weiland, J. L., Hauser, M. G., & Kelsall, T. 

1997, ApJ, 480, 173 
Sparks, W. B., Macchetto, F., & Golombek, D. 1989, ApJ, 345 153 
Temi, P., Brighenti, F., Mathews, W. G., k Bregman, J. D. 2004, ApJS, 151, 237 
Temi, P., Brighenti, F., & Mathews, W. G. 2007, ApJ, 660, 1215 
Temi, P., Brighenti, F., & Mathews, W. G. 2008, ApJ, 672, 244 

Tielens, A. G. G. M., McKee, C. F., Seab, C. G., & Hollenbach, D. J. 1994, ApJ, 431, 321 
Tielens, A. G. G. M. 2008, ARA&A, 46, 289 

Tonry, J. L., Dressier, A., Blakeslee, J. P., Ajhar, E. A., Fletcher, A., Luppino, G. A., Metzger, M. 

R., & Moore, C. B. 2001, ApJ, 546, 681 
Trinchieri, G., Pellegrini, S., Welter, A., Fabbiano, G., & Fiore, F. A&A, 364, 53 
Tsai, J. C. & Mathews, W. G. 1995, ApJ, 448, 84 
Verdoes Kleijn, G. A. & de Zeeuw, P. T. 2005, A&A, 435, 43 

Verdugo, E., Yamamura, I., & Pearson, C. P. 2007, AKARI FIS Data User Manual Version 1.2 

Werner, M. W., et al. 2004, ApJS, 154, 1 

Wiklind, T., Combes, F., & Henkel, C. 1995, A&A, 297, 643 

Wiklind, T. & Henkel, C. 1995, A&A, 297, L71 

Willner, S. P., Elvis, M., Fabbiano, G., Lawrence, A., & Ward, M. J. 1985, ApJ, 299, 443 
Xilouris, E. M., Madden, S. C., Galliano, F., Vigroux, L., & Sauvage, M. 2004, A&A, 416, 41 



25 



